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Abstract 

The key problem of the theory of exterior differential systems (EDS) is 
to decide whether or not a system is in involution. The special case of EDSs 
generated by one-forms (Pfaffian systems) can be adequately illustrated by 
a 2-dimensional example. 

In 4 dimensions two such problems arise in a natural way, namely, the 
Riemann-Lanczos and the Weyl-Lanczos problems. It is known from the 
work of Bampi and Caviglia that the Weyl-Lanczos problem is always in in- 
volution in both 4 and 5 dimensions but that the Riemann-Lanczos problem 
fails to be in involution even for 4 dimensions. However, singular solutions 
of it can be found. 

We give examples of singular solutions for the Godel, Kasner and Debever- 
Hubaut spacetimes. It is even possible that the singular solution can inherit 
the spacetime symmetries as in the Debever-Hubaut case. We comment on 
the Riemann-Lanczos problem in 5 dimensions which is neither in involution 
nor does it admit a 5-dimensional involution of Vessiot vector fields in the 
generic case. 
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I. Introduction 



The problem of generating the spacetime Weyl conformal curvature tensor 
C a bcd from a tensor potential is called the Weyl-Lanczos problem and the 
analogous problem for the Riemann curvature tensor the Riemann-Lanczos 
problem. 

The Lanczos tensor potential admits the following index symmetries 

L[ab]c = L a b c , (1) 

where a, b,c,s = 0, 1, 2, 3 and 

L[ a bc] = . (2) 

Apart from these, we may impose two gauge conditions: the differential 
gauge condition 

L a b s - S = , (3) 

where ";" indicates covariant differentiation, and the algebraic gauge or trace 
free condition 

L a s s = . (4) 



Lanczos discovered the Weyl-Lanczos equations [3C], where he intro- 
duced a Lagrangian based on the double dual of the Riemann tensor R a bcd- 
The Lanczos tensor L a i, c arose as a Lagrange multiplier for this Lagrangian. 
Lanczos found an expression for the Weyl tensor in terms of certain Lagrange 
multipliers L a fcc, namely, 

Cabcd = L a bc\d ~ L a bd;c + L c da;b ~ L c db;a + 9bcL( a d) + 9adL(bc) ~ 9bdL( a c) 
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—9acL(bd) + -L mS m;s(9ac9bd ~ 9ad9bc) , (5) 

where L a d = L a S d- s ~ L a s S \d- These are the Weyl-Lanczos equations. They 
can also be expressed in the form 

Cabcd = L[ab][c;d] + L[ cd ][ a - b ] - £[ ab ][ c;(i ] - £* cd ] [ a; fe] (6) 

as in 37]. Further, if we impose (||) and (§) above, we can simplify (||) 
considerably to 

Cabcd = L a bc\d ~ L a bd;c~\~ L c da;b ~ L c db\a ~ 9bcL a d-s ~ 9adLl cs +9bdL a c-s +5ac-^6d- s • 

(7) 

Many solutions to (^) are known and solutions for vacuum spacetimes can 
be found in [13]. 

We can also attempt to express the Riemann curvature tensor in terms 
of a comparable tensor potential L a bc which leads to the Riemann-Lanczos 
problem for spacetimes and which has been discussed at length by Bampi 
and Caviglia ||, |J. Udeschini Brinis Q had wanted to describe the space- 
time Riemann tensor in terms of a Lanczos tensor L a b c an d proposed the 
Riemann-Lanczos relations 

Rabcd = L a bc;d ~ L a bd;c + L c d a - ~ L c db- a ■ (8) 

But the difficulties with the relations @ were pointed out in two papers by 
Bampi and Caviglia ||, 0] where they proved existence theorems for solutions 
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of (||) and (||). Whereas the Weyl-Lanczos problem is always in involution, 
the Riemann-Lanczos problem is not and only singular solutions of it can 
occur if the problem is not modified. Bampi and Caviglia showed that: 

i) The Weyl-Lanczos problem (||) or (j^) has non-singular solutions for n = 
4, 5, where n is the dimension of the space-time manifold M ; 

ii) For n = 4 the Riemann-Lanczos problem (|8|) has no non-singular solutions 
but it does have "singular" solutions which means that the Cartan characters 
do not adopt their maximal values; 

iii) The differential gauge condition (^) has no effect on the existence or 
non-existence of solutions of either (||) or (||). 

Bampi and Caviglia Q also suggested a prolongation of the Riemann- 
Lanczos equations to a second-order system which makes them a system in 
involution. 

In the Riemann-Lanczos problem, we meet equation (||) and possibly (||) 
but not (Q), which leaves us with 20 independent components for the L a bc 
in 4 dimensions. We always assume that the cyclic conditions @ hold but 
not the trace-free conditions (||). If equations (Q) were to hold, we would 
have 

R = 4L nk k . n = -AL nk n . k = , 

which would lead to inconsistencies. Because we are only going to talk 
about the Riemann-Lanczos problem in this paper, we shallchange notation 
from L a h c to L a h c from now on. We shallalso write the Riemann-Lanczos 
equations in solved form as 

fabcd := Rabcd ~ ^abc\d + ^abd\c ~ L c da;b + ^cdb\a ■ (9) 



The paper by Massa and Pagani |32j concerning a modification of the Riemann- 
Lanczos problem used a totally different approach to the above work of 
Bampi and Caviglia. Accordingly, we do not consider |33] here. As refer- 



ence [15] is based on [32] it is also not applicable here. 



II. Exterior Differential Systems and Cartan Char- 
acters 

It is necessary to introduce some notions from the theory of exterior differ- 
ential systems (EDS) in order to understand the structure of the Riemann- 
Lanczos problem. We denote a formal iV-dimensional manifold by A4 of 
which the spacetime manifold M is an n-dimensional submanifold. Then, 
for a collection of differential forms we define an exterior differential system 
(EDS). 

Definition 1 An exterior differential system EDS £ on M consists of 

1 — forms 

2 — forms 

p — forms afj (1 < i p < k p ) , 
and sometimes 0-forms on Ai. 

Any EDS E is a subset of the set of all differential forms on M. There is 
a further set of differential forms on Ai formed by a subset of £ called its 



„(1) 



ai' (1 < k < h) , 

(2) 
Oi- 

12 



(l<h< k 2 ) 
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associated ideal 2(T,) 0, The associated ideal Z(E) is a differential 
ideal X, that is, dZ C X. This means that such an ideal is closed under 
exterior differentiation. If the exterior differential system (EDS) £ itself 
is closed under exterior differentiation, we say that £ is a closed exterior 
differential system. Next, we wish to look at Pfaffian systems and Vessiot 
vector fields as they are of great importance to our applications later on. 



A. Pfaffian Systems and Vessiot Vector Fields 

A Pfaffian system P is a special EDS containing only one-forms 9 a and 
0- forms. We denote the rank of V by s which is given by the number s of 
independent one-forms in V . Consider now the collection of all vector fields 
on J\4, X(M). There is a subset of X{J\A) given by all those vector fields 
which annihilate the Pfaffians in P 

T>:={Xe X(M)\9 a (X)=0, a = l, ••-,«} 

so that the number of independent vector fields in T> is N — s. We call T> 
the dual system to P. For any Pfaffian system V, we can also look at its 



derived system V' for which we can write [41] 



P' := {0 a € V\d6 a (X, Y) = 0, MX, Y G V, a = l,---,s}, 

where V CP always holds. 

Pfaffian systems are classified according to the ease with which they can 
be integrated. The most familiar Pfaffian Systems are also the simplest, 
namely, the complete Pfaffian systems to which the celebrated Frobe- 
nius theorem applies. We say that a Pfaffian system P is complete if 
P = P' holds. A necessary and sufficient condition for a complete Pfaffian 
system is now given by the Frobenius theorem. 

Theorem 1 Frobenius Theorem 

Let the s independent 1-forms Q^,...,Q S generate a closed ideal I (V), then, 
we can find coordinates in a neighbourhood of any point x l , ...,x such that 
all forms are generated by the s coordinate differentials dx N ~ s+l , • • • , dx . 



A good account of Pfaffian systems can be found in [21, g] or in a more 
theoretical approach in [ 20 1 . In general, Pfaffian systems are not complete 
but each Pfaffian system can be enlarged by adding further 1-forms until the 
enlarged system is complete or the enlarged system becomes inconsistent. 
The minimal enlarged Pfaffian system which is complete is known as the 
associated system A(P) of P and its dimension is called the class c of P. 
Sometimes, A(P) is also called the Cartan system C(P) p8[| . 

A Pfaffian system P can dually be characterised using Vessiot vector 



field systems as discussed in ]41j]. This approach was developed in Vessiot 



[44] and can be found in [18]. First, we take a vector field system V, which 
we here choose to be the dual system T> of a Pfaffian system P so that V = T> 
is given by N — s vector fields Y l and we define in turn its first and second 
derived systems as 



V := V+ [V, V] 
V" := V + [V ,V] 



(10) 
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and so on [41|. For the derived system, the inclusion T> C T>' is always valid 
and (T 3 ') = T>' for a Pfaffian system, where T> = . We can also formulate 
the notion of completeness in terms of vector field systems and versions of 
the Frobenius theorem in this language are numerous. We therefore define 
the space V C T> as 

T> := {Y E V\d6 a (X, Y) = VX 6 V, a = 1, • • • , s} . 

Here, we state the dual version of the Frobenius theorem as 

Theorem 2 Complete Pfaffian Systems 

Given any two vector fields X, Y G V, then V is a complete Pfaffian 

system means that the the commutator vector field [X, Y] is also a vector 
field in T>. 

A Vessiot vector field system T> is complete means then that [T>,T>] C T> for 
all Y € T> which is equivalent to T> = T>. 

When a vector field system T> or Pfaffian system V fails to be complete, 
a slightly weaker condition may hold on some subsystem of T>. Such a 
possibility is that a vector field system is in involution, which is defined as 
follows: 

Definition 2 Involutory Subsystem of a Vector Field System 
An involutory subsystem T of a vector field system T> is a subsystem T of 
T> such that [T, T] C T>, which means that T is closed relative to T> but not 
relative to T itself. 

This means that if we are given two vector fields X,Y € T, it is then valid 
for all 6 a G V 

d6 a (X,Y) = 9 a (X)Y -9 a (Y)X -^9 a ([X,Y]) 

= -\e a ([x,Y]) 

+ 0. (11) 

This is because we have 9 a (X) = and 9 a (Y) = holding identically but 
not \9 a {[X,Y)) = since [X, Y] need not lie inside T. However, if we allow 
for the vector fields in T> — T, then 9 a ([X, Y]) can vanish i.e. 

d6 a (X,Y) = 0(modJ(V)). (12) 

Note that when we take T = T>, we get a complete vector field system but 



when T is strictly contained in T> it can fail to be complete [41]. 

Further, the dual space of C(V) which is usually denoted by C{T>) = 
C(V) 1 ' is the space of all Cauchy characteristic vector fields, where such a 
vector field is defined as 

Definition 3 Cauchy Characteristic Vector Field 

Y is a Cauchy characteristic vector field of a vector field ideal I(V) means 
that 

Y\6 a = 
Y\dO a = A 7 # 7 = {mod V) 
for all one-forms 9 a in V , where A 7 are scalar multipliers for each 6 a . 
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Both C(T>) and C(V) are complete Pfaffian systems and it is C(T>) = T>. 
For their dimensions it has to hold dim(C(P)) + dim(C('P)) = N. Cauchy 
characteristic vector fields can also be characterised as || 

Theorem 3 Given a Pfaffian system V and its associated ideal T{V) and let 
T> be the dual space ofV . Y £ T> is Cauchy means that [X, Y] 6 D VX € T>. 

For EDS and Pfaffian systems in particular, it can be of interest to deter- 
mine its symmetries and especially their infinitesimal generators. In the 
same way as Killing vector fields are symmetries of the metric tensor and 
hence of Einstein's equations, any EDS or system of PDEs can adopt sym- 
metries. When an EDS adopts symmetries it means that there exist vector 
fields Y whose Lie derivatives applied to each form in the EDS annihilates 
them or those new forms are contained in the EDS itself. These infinitesimal 
symmetry generators Y are called isovectors and defined as 

Definition 4 Isovectors 

Y is an isovector of an EDS £ means that £yS C £ and dually Y is an 
isovector of a vector field system V means that £yT> C T>. 

Cauchy characteristics form a special sub-algebra of all isovectors of a system || . 
Isovectors Y do not necessarily have to be such that Y £ T>. Often there 
exist isovectors Y such that Y ^ T>. But for those Y with Y £ D, it is 
known that || 

Theorem 4 Cauchy Characteristics 

Given a vector field ideal T(V) . Then, Y € T> is an isovector means that Y 
is a Cauchy characteristic vector field. 

Note that if V is complete, then all its vector fields Y € T){= T>) are isovec- 
tors (= Cauchy characteristics). The formula of H. Cartan 

£ Y a = Y\da + d(Y\a) (13) 

is used when we determine the isovectors of p- forms a. See R, [|, [L4|, [5], 
|35f for discussions of symmetries and isovectors. An algebraic computing 
package assisting with symmetry determination is the REDUCE package 
DIMSYM by Sherring 

We illustrate the above theory with the example of a coordinate trans- 
formation to Euclidean coordinates in two dimensions. For the line element 
we therefore have 

ds 2 = y 2 dx 2 + dy 2 = du 2 + dU 2 , 
where u, U are the Euclidean coordinates. For this to hold we must have 

y 2 = p 2 + p 2 , 

= pq + PQ, 

1 = q 2 + Q 2 , (14) 

where we introduced the Monge notation 

_ du _ du _ d 2 u _ d 2 u _ d 2 u 
^ dx ' ^ dy ' dx 2 ' dxdy ' dy 2 ' 

P .= ^L n -=^L b - d2u «? - d2jJ r-=^E 

dx dy dx z dxdy dy z 
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These local coordinates form a jet bundle l 7 2 (M 2 ,R 2 ) with N = 14 for- 
mal dimensions, where we have n = 2 independent and m = 2 dependent 
variables. After differentiating equations (|l4] ) with respect to x and y and 
rearranging them, we obtain a Pfaffian system V with s = 12 independent 
1-forms which can locally be expressed as 

9 1 = du — pdx — qdy 
dU - Pdx - Qdy 
dp — rdx — sdy 
dP — Rdx — Sdy 
dq — sdx — tdy 
dQ — Sdx — Tdy 
pdr + rdp + PdR + RdP 
qdr + rdq + QdR + RdQ + dy 
pds + sdp + PdS + SdP - dy 
qds + sdq + QdS + SdQ 
pdt + tdp + PdT + TdP 

qdt + tdq + QdT + TdQ . (15) 

As the rank of T> is N — s = 14 — 12 = 2, the system (|l5| ) can dually be 
characterised by two Vessiot vector fields V 1 ,V 2 generating T>, where 

V* = V^ + V^ + ( P V: + qV^ + (PV: + QV^ + (rV: 



9 9 

,10 
)12 



d d ■ d ■ d 

(SV l + TV 1 )— + V — + VI— + V — + V 



+ ^ + ^' < = 1 ' 2 (16) 

and the coefficients VJ?, V^, V*, Vg, V t l , Vp are determined through the con- 
dition that the six 1-forms 6 7 to 12 are annihilated by V 1 ,V 2 . Therefore, 
we obtain 

K = h&(rs + RS)-r 2 -R 2 )Vi + (^(rt + RT + l)-rs-RS)Vj], 
a Q Q y 

V l R = -^-[{rs + RS + -(^(rs + RS-r 2 -R 2 ))V* + (rt + RT+l 
Q a Q 

+~(£(rf + RT + 1) - rs - RS)Vj] , 

V* = h^(s 2 + S 2 )-rs-RS)Vl + £(st + ST)-s 2 -S 2 + l)V*], 
a Q Q y 

V l s = -h(s 2 + S 2 +«£(s 2 + S 2 )-rs-RS))Vj: + (st + ST 
Q a Q 

+^(st + ST)-s 2 -S 2 + l)Vj], 

V * = \[^ st + ST )- rt - RT )n + {^{t 2 + T 2 )-st-ST)V^], 

Vi = -h( s t + ST + ^(st + ST)-rt-RT))V^ 
Q a Q 



+ (f + T 2 + -&{t 2 + T 2 ) -st- ST))]Vj , (17) 
a Q 
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where a is given by a = pQ — Pq. Equations ( |17| ) tell us that all components 
of V , V 2 are determined by V x , Vy which can be assigned arbitrarily. In this 



way all 2-forms arising from (15) and applied to V , V 2 vanish identically. 
This means that the system ( |T5| ) is complete and its derived system consists 
of the two Vessiot vector fields V 1 , V 2 spanning T> so that T> = T> = T>' . It 
also means that C(V) = T> so that c = 2 and therefore, both V 1 and V 2 
are Cauchy charcateristics. At the same time, they form an involution of 
maximal dimension g = 2. Note that the V 1 , V 2 are also isovectors but that 
other isovectors which are not in T> such as Y l = C l u + Cy-J^y n= 1,2 do 
occur. 

B. Integral Elements, Cartan Characters and Genus 



All definitions in this section are based on [33] which relies on Cartan's origi- 
nal work whereas the work in |^7| varies slightly from [ 10 1 in some definitions. 
When an EDS £ or a Pfaffian system V is given, we try to find manifolds on 
which all the differential forms of £ or V are annihilated. Firstly, we shall- 
construct tangent spaces to such manifolds, where the dimension of each of 
them can be deduced from a sequence of non-negative integers called the 
Cartan characters. 

Assume, we are given a vector space M N of which we want to find all 
possible distinct p-dimensional subspaces. The set of such subspaces forms a 
p(N — p)— dimensional manifold called the Grassmann manifold G(N,p). 
We can characterise G(N,p) as a quotient space of orthogonal groups as 
follows: 

G{N,p) ~0(N,p)/0(N)xO(p) . 

Alternatively, one can write each p-dimensional vector space as a p-vector 
v = Ae 1 A • • • A e p , where {e 1 , • • • , e p } is its basis. Using this definition we 
obtain 

G(N,p)~{(v l ,---,v p }/~, 

where (v 1 A • • • A v p ) ~ (w 1 A • • • A w p ) 43> (to 1 , w p ) = A(V, • • • , v p ) 

meaning that all multiples of a p-vector are in the same equivalence class. 
The Pliicker coordinate of a p-vector v = ■ ■ • ai p e n A • • • A e %v is 
given by the equivalence class [(a^ • • • a^,)] of which each equivalence class 
represents a distinct p-plane. 

Because we are not only looking at a single G(N,p) but at a whole set 
for a given manifold M, we define the Grassmann bundle as 

G P {M):= |J G(N,p)\ x , 

where G(N,p)\ x is the Grassmannian manifold at a point x G Ai. We 
shalluse these definitions later on to help to describe integral manifolds but 
first we introduce integral elements as 

Definition 5 Integral Element 

A p-dimensional subspace (E p ) x of the tangent space T X (A4) at a point x 
on a N -dimensional manifold Ai is called an integral element of an EDS E 
if a l j(E p ) x = at x for any form a l j in S p , 1 < i < p, 1 < j < k p , which 
means that all differential forms o/S are annihilated on (E p ) x at x £ A4. 

Each p-dimensional integral element {E p ) x is then a particular element of 
G(N,p)\ x . The set of all p-dimensional integral elements (E p ) x for all x G Ai 
is then denoted by V P (S), where V P (S) C G p (Ai). 
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Once we have found such a p-dimensional integral element (E p ) x C T X A4 
spanned by V 1 , ■ ■ • , V p , we are looking for vectors V£ +l € T X M. in such a 
way that the space generated by (E p ) x and V x +1 form a (p + l)-dimensional 
integral element. The conditions on such a tangent vector V^ +1 at the point 
x £ M. are: 

= 0, 1 < ii < fei , 

= 0, l<i2<k 2 ,l<j 1 <p, 



^(V 1 ^,---^^ 1 ) = 0, l<i p+1 <k p+1 . (18) 

The vectors V^" 1 " 1 that satisfy the above polar system H((E P ) X ) of linear 
equations generate the polar space H((E P ) X ) ± of (E p ) x . Depending on 
the ranks of the polar systems generated, we can divide integral elements 
into 3 classes which we call regular, ordinary and singular. We look at 
the subsystem a}(V E 1 ) = ■ ■ ■ = a\ {V x ) = of (18) where sq{x) is its rank. 
From this, we get the integer 



s := max {x( z M }(r(H((E ) x )) = s (x)) 

called the zeroth Cartan character. We define a regular point iEjM 
to be a point where sq(x) = Sq. Then, a 1-dimensional integral element 
{E v ) x is an ordinary integral element if x is a regular point. Let the polar 
system of (E l ) x in T X M have rank r(H((E l ) x )) = s\{x) + sq. From this, 
we define 

si := max {xGM y2 £TxM} r(H((E 1 ) x )) - s , 

which is called the first Cartan character. If, for x G A4, s\(x) = s± 
holds, then the integral element (E 1 ) x is called regular. A 2-dimensional 
integral element (E 2 ) x then is called ordinary if it contains at least one 
regular 1-dimensional integral element. We define inductively 

Definition 6 Cartan Characters 

The p th Cartan character is inductively defined as: 

s p = max {x€Myp+ i €TxM} r(H((E p ) x )) - Y^Zl «»• 

Using the notion of Cartan characters we can now precisely define 

Definition 7 Regular, Ordinary and Singular Integral Elements 
A p- dimensional integral element (E p ) x is called ordinary if it contains at 
least one (p-1)- dimensional regular integral element. A p-dimensional inte- 
gral element is called regular if its polar system H{{E P ) X ) has maximal rank 
r(H((E p ) x )) = s p (x) = s p . A p-dimensional integral element is singular if 
it is neither regular nor ordinary. 

Note that a sequence of integral elements {E°) x C {E l ) x C ••• C {E p ) x 
at a point x € Ai is called a regular chain of integral elements if all its 
(E k ) x , 1 < k < p — 1, are regular and (E p ) x is at least ordinary. The 
maximal dimension an ordinary integral element can adopt is called the 
genus g of S or V . 

Integral elements are the tangent planes to manifolds which are solution 
manifolds for a given EDS E and Pfaffian systems V and defined as 
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Definition 8 Integral Manifolds 

An integral manifold of an EDS £ on A4 is a p- dimensional submanifold 
M of M such that each k-dimensional vector subspace (E k ) x C T x J\f for 
1 < k < p annihilates all the k- forms in the EDS S. 

Next, we wish to illustrate this theory using example (|l5| ) again. There, a 
1-dimensional integral element (E 1 ) x can for instance be created by means 
of setting V x := 1, V y l := and computing the other components according 
to (|j~6|) and ([H]). Its polar space is then spanned by H{(E l ) x ) = {V 2 } for 
any V 2 satisfying ( |l6| ) and ( |T7| ) other than V 1 itself. Then, (E 2 ) x is given 
by (E 2 ) x = {V 1 , V 2 } which is the maximal dimension an integral element 
can adopt for this example, where we can for instance choose V 2 := and 
V 2 := 1. The polar matrix of the polar space H(E°) X ) is given by 



dx 


dy 


du 


dU 


dp 


dP 


dq 


dQ 


dr 


dR 


ds 
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dt 
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-Q 
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The rank r(H(E°) x ) of the polar matrix is equal to the rank s of our Pfaffian 
system (|15|) and therefore we obtain so = s = 12. The two Vessiot vector 
fields V 1 which form a 2-dimensional involution for (jl^) also constitute a 
regular chain of integral elements where (E l ) x is spanned by either of the 
V 1 and (E 2 ) x by both of them. 

Next, we wish to determine the integral manifold for ( |l5| ) on which ( |14| ) 
holds as well. Prom the 6 equations obtained from differentiating ([lj) with 
respect to x and y, we see that t = = and T = = so that the 
solution must be linear in y 

u = yf{x)+g{x) 

U = yF{x) + G{x) (19) 

with f(x),g(x),F(x),G(x) arbitrary functions. From this we find by very 
elementary calculations that the general solution is given by 



u 


= y cos (x + e) + a , 


U = 


y sin (x + e) + b 


p 


= — y sin (x + e) , 


p = 


y cos (x + e) , 


(I 


= cos (x + e) , 


Q = 


sin {x + e) , 


r 


= -y cos (x + e) , 


R = 


y cos (x + e) , 


s 


= — sin [x + e) , 


S = 


cos (x + e) , 


t 


= 0, 


T = 


0. 
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where {a, b, e} are arbitrary constants corresponding to the translational and 
rotational degrees of freedom. This also gives a local parameterisation of 
the 2-dimensional integral manifold which corresponds to this solution. 



C. Independence Condition, Reduced Characters and Invo- 
lution 

Until now, all variables of a given EDS £ were treated equally. But if we wish 
to look for specific integral manifolds transversal to some given submanifold, 
we must specify this fact for our EDS. This is achieved by introducing an 
independence condition 17 = uj 1 A • ■ ■ A uj n 7^ 0, which is sometimes also 
called the volume element, where the lo 1 are 1-forms characterising such a 
submanifold. If we select local coordinates, where we shalluse 

brackets to indicate powers of x l such as in (x t ) n , then 17 is given by 17 = 
dx 1 A • • • A dx n . 

All integral elements on which 17 does not vanish are called admissible 
integral elements according to [45] and using this concept we can define 



Definition 9 Involutive Systems 

An EDS with independence condition (£, 17) is in involution with respect 
to 17 at a point x £ Ai means that there exists an admissible ordinary 
integral element of dimension n of (E, 17) at x. 

For practical reasons new integers called the reduced Cartan charac- 
ters are introduced so as to decide whether a given EDS E is in involution 
or not with respect to an independence condition 17. In order to deter- 
mine them, we need to introduce the reduced polar systems H red ((E p ) x ) 
which are defined as the polar systems H red ((E p ) x ), where all terms in- 
volving any of the u , ■ ■ ■ ,u> n are suppressed. We define the reduced rank 
s' (x) := r(H red ((E°) x )) and s' := rnax^^^s'^x)) so that we can define 
the reduced characters inductively. 

Definition 10 Reduced Cartan Characters 

The reduced Cartan characters s' p are inductively defined as 



s' p := rnax( [xeMyr > eTxM} r(H^ d ((E p ) x ))) -J^4 

i=0 



The coincidence of the set of characters with that of the set of reduced 
characters is a necessary condition for a system to be in involution. In 
practice, we can use Cartan's test for involution based on a comparison 
of polar elements with the reduced polar elements. But first, we must enlarge 
our account of Pfaffian systems with the use of more general concepts from 
the literature. Given a Pfaffian system V consisting of s 1-forms 8 a with 
independence condition 17 = lo 1 A ■ ■ ■ A uj n 7^ 0, we denote by 7r A all the 
extra forms such that (9 a , u i , tt x ), 1 < a < s, form a coframe on our formal 
iV-dimensional manifold Ai. Once we have chosen such forms 7r A , we can 
write each d9 a as 

d6 a = A^tt x A J + -B^J A J + -C^V A tt k (mod 1{V)) . (20) 



In equations (p0|), the A^ form the tableau matrix and the Bf- are called 
the torsion terms 1 . The Pfaffian system is called quasi-linear if all the 
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C? = 0. The following changes of coframe preserve linearity 



We say that the torsion can be absorbed if a suitable transformation <E> 

$ : 7T A -> 7T A + P^W* 



= B?. + - A a x( p) (22) 



with 

can be found such that 5?- = 0. When Bf- ^ in every coframe (6 a , uo e , tt x ), 
the system possesses integrability conditions which will prevent our given 
Pfaffian system V from being in involution. If a transformation <I> above 
exists such that all above Bf- = 0, then there are no integrability conditions 
to stop the system being involutive. 

We can rewrite ( |20| ) in the dual language of Vessiot vector fields. If 
we choose a frame dual to the coframe (u> 1 , 9 a , 7r A ), which we denote by 
(Ei,F a , G\) such that 





= 




■ ,n) 






= d f3 ) 


(a, /3 = 1, • 






tt x (G k ) 




(A,« = 1, • 


■■N-n-s) 


(23) 



two Vessiot vector fields 1,7 £ P can then be recast as 

X = X^Ei + XfyF a + X^Gx , 

Y = %) E i + Y (6) F P + Y (l) G K- ( 24 ) 
When we calculate d8 a (X, Y), we obtain that 

d9 a (X, Y) = A\\(X^Yfa - X^Y^) + BfjX^Y^ + C{* K X^Yfo 
= mod (1(V)) , (25) 

where (p3) is a system of equations for X,Y 6D based on the dual approach 
to (120). 



From the tableau matrix A?-, we can determine the reduced Cartan 



characters directly [j45| , where for the ranks r(^4° 1 ) := s[ and so on. 
There are many versions of Cartan's original test for involution [1C] some of 
which one can find in || or in |^], where a version for Pfaffian systems is 
discussed. We state it as: 

Theorem 5 Cartan's Test for Involution 

A Pfaffian system (V, fi) is in involution means that sq = s' , s\ = s' 1; • • • , s p = 
s' p and a coframe transformation can be found such that all its torsion terms 
Bfj vanish identically. 

The involutive genus g determines the maximal dimension a regular in- 
tegral manifold can adopt and is given by g = N — ^7=0 s i- Note that 
the notion of an involutory subsystem of a vector field systems is not the 
same as the one used in discussing differential forms. In order to obtain 
the equivalent of this, we must add an independence condition and define 
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what it means for a vector field system to be in involution with respect to 



x 



The precise definition can be found in [41]. We can also use 



Groebner bases in developing criteria for involutivity. This is carried out in 
many places for instance in a paper by Mansfield [31]. 



In our example (15), we choose w 1 := dx and w 2 := dy so that Q = 
dx A dy whereas the u and U constitute the dependent variables. Then, the 
reduced Cartan characters happen to coincide with their full counterparts so 
that s = so = s' = 12 while all higher reduced characters vanish. Because 
the 12 6 a in (|i~5| ) together with uj 1 = dx and uj 2 = dy form a complete 
coframe of dimension N = 14 already, we do not need to add any ir x and 
the tableau matrices vanish identically with all components = 0. This 
proves that s^ = 0, \/i > 0. The two torsion terms B\ 2 and B 2 2 vanish 
identically. Next, we must compute the torsion terms which do not vanish 
identically and they are 

Bf 2 = B\ 2 = B\ 2 = Bf 2 = rt + RT + 1 - s 2 - S 2 . 

But for our system transforming polar into Cartesian coordinates we have 
t = T = and s 2 + S 2 = 1 so that all the torsion terms vanish identically 
and the system is in involution, and, the involutive genus g is given by 
g = 14 - 12 = 2. 

D. Existence Theorems and Prolongation 

An important question is whether and when integral manifolds exist and 
whether, given some initial data, they are unique. The Cartan - Kahler 
theorems ||, |45| give some answers to these questions for real-analytic 
systems. They reduce to the Cauchy-Kovalevskaya theorem when the 
number of equations equals the number of unknowns which is given by m. 



A version for first-order systems is given in [23] which is sufficient because 
every system of PDEs can be transformed into a first-order system of PDEs. 

The first Cartan-Kahler theorem specifies under which conditions a unique 
integral manifold of dimension p can be constructed from a (p— l)-dimensional 
one. For the second Cartan-Kahler existence theorem, we need our EDS E 
to be given in normal form. Then, the second theorem simply states that 
under certain conditions a whole chain of regular integral manifolds of in- 
creasing dimensions exists (see Appendix A). 

If there is no general solution for a given EDS (E, f2), then there could 
be identities or even integrability conditions so that we must prolong the 
system by adding in these conditions and by enlarging the space of our formal 
jet coordinates to higher order. Good instructive examples of prolongation 
methods were given in [|ll], 17] and an example of an infinite dimensional 



prolongation of vector fields can be found in [|19[. Next, we shalldiscuss 
Cartan's classical approach for an EDS briefly. But if we are dealing with a 
system of PDEs, then prolongation means just adding the partial derivatives 
of order q + 1 as new jet-coordinates and supplementing our original system 



of PDEs with all partial derivatives of the equations in the system [35], [38]. 



Cartan's classical procedure [1Q], |J or [45] starts with an EDS with 
independence condition (E, O) given as 

E = (Fx,.. ., F ko ), k - forms 
Ei = (af\ ...,a^), h 1 - forms 
E 2 = (af\ ...,af 2 ) ), k 2 2 - forms 
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(26) 



and so on. We now complete {a^\oj l ) to a coframe {a^\oj % ,ir x ) on M.. 
We then choose an admissible integral element (E n ) x from the set of all 
admissible integral elements over A4, which we denote by V(E, f2) C V n (E) , 
on which 

qJ' = ti u iJ = and vr A = t,V 

has to hold for some (t^ which are our new coordinates on the fiber 
G P {M) X = G(N,p) , where T X {M) = W. Based on this we can define the 
first prolongation of (S, f2) as 

Definition 11 First Prolongation 
The first prolongation E^ of (E, fi) is formed by 
i) the original EDS (E, Vt) restricted to V(E, f2) and 
iij r/ie closed system on V(E,0) given by 



and 



7T = 7T — Ej (J , 



We essentially added an extra number of contact conditions and their ex- 
terior derivatives to the above EDS (E,0). Higher-order prolongations are 
defined accordingly and for a detailed account on prolongations see @ . The 
key result for prolongations in the real analytic case [29] is given by the 
Cartan-Kuranishi theorem 

Theorem 6 Cartan-Kuranishi Theorem 

An EDS (E,Q) which is not in involution and not inconsistent becomes 
either involutive or inconsistent after a finite number of prolongations. 

An upper limit for the number of prolongations for involutivity or inconsis- 
tency is given by a number q which can be determined inductively (c.f. |36[| , 
J38[ | or pG|]). In the next sections, we shallapply some of the above theory 
to the Riemann-Lanczos problems in 4 and 5 dimensions. 



Concluding Remarks: We should point out that there is a quantity simi- 
lar to the Cartan characters also introduced by Cartan [10] called the degre 
d'arbitraire. Seiler derived a relation between this degree of arbitrariness 
and the Cartan characters for nonlinear, overdetermined higher-order sys- 
tems of PDEs [39]. Einstein introduced the strength of a system of equa- 
tions in his last unified field theory of electromagnetism and gravity in order 
to measure how strongly such a system restricts the solution. A discussion 
between him and Cartan on this subject is presented in . Sue finally com- 
pared Cartan's degree of arbitrariness and Einstein's strength in [42, Q for 
linear systems of equations. 



III. The Riemann-Lanczos Problem as an EDS and 
Singular Solutions 

In the Bampi-Caviglia papers ||, |j the Riemann-Lanczos problem was writ- 
ten as an EDS and found not to be in involution in 4 dimensions as mentioned 
earlier. In their second paper [fl they introduced a prolongation to create an 
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involutive system. Now, we are going to form the Pfaffian system resulting 
from the Riemann-Lanczos equations (||). The Pfaffian system we are go- 
ing to introduce slightly differs from that in || because we incorporate the 
cyclic conditions (|2|) as well. The system consists of 20 exterior derivatives of 
the Riemann-Lanczos equations together with the 20 contact conditions and 
their exterior derivatives and we can omit the differential gauge conditions 
which do not change our results qualitatively. Using the exterior differen- 
tials of equations @ in solved form together with the contact conditions, 
we obtain the system 

dfabcd = (Rabcd,e + CXabcde)dx e — dP a bcd + dP a bdc ~ dP c dab + dP c dba 
+^ad(dL n bc + dL nc b) — T^ c (dL n b d + dL n db) 

~^^bc(dL na d + dL n da) — ^bd(dL nac + dL nca ) 

K a bc — dL a b c Pabcedx 

dK abc = dx e AdP abce , (27) 
where a is given by 

Ot-abcde = ^ad,e(Pnbc + L nc b) — T^ ce (L n b d + L n db) + T bce (L nac i + L n d a ) 
~ L bd,e\^nac 

+ L 

■nca ) ■ (28) 

The (V, L abc , P a bcd) are local jet coordinates on ^(M 4 , M 20 ) with P a b c d = 
d g^d c when projected onto our spacetime manifold, and the K a b c are the 20 
contact conditions. A Vessiot vector field V € V is locally given as 

Q Q Q 

V = V e — + V e P abce —— + V abcd —^, (29) 

dx e dL abc dPabcd 

where 20 of the 80 components V a b c d are defined through the requirement 
that 

dfifcdiV) = • (30) 
This leads to the 20 independent equations 

V{abcd} = y £ [Rabcd,e + Oi a bcde +^ad(Pnbce + Pncbe) —^ac(Pnbde + Pndbe) 
~^~^bc(Pnade Pndae) F bd(Pnace ~\~ Pncae)] > (31) 

where {abed} denotes the Riemann-type symmetrisation performed over the 
indices abed. The 60 remaining components Vabcd — V{ a b c d] can be chosen 
arbitrarily. As shown in || a rank deficiency between the Sj and the s- 
occurs leading to 6 identities and the system (27) fails to be involutive. A 



more detailed explanation on how those identities can be created as linear 
combinations of rows of polar space matrices can be found in [[22j| . 

The reduced Cartan characters can be computed using the REDUCE 



computer code in [22] based on the EDS package by Hartley [24] and slight 
modifications of it. For a number of spacetimes such as conformally flat, Kas- 
ner, Godel, some Debever types and plane wave spacetimes we verified that 
the reduced Cartan characters (s' , s[, s' 2 , s' 3 , s' 4 ) were (40,20,19,15,6) when 
the differential gauge condition (|3|) was not included and (40,20,19,15,0) 
when it was included. In both cases, the systems were clearly not in involu- 
tion because some of the reduced characters differ from the Cartan charac- 
ters - a fact reflected by results given by the computer code. 

Even though the Riemann-Lanczos problem in 4 dimensions is not in 
involution, singular solutions do exist. Singular solutions occur when some 
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of the integral elements, which are tangent to an integral manifold, are 
singular leading to singular integral manifolds, which means that the Cartan 
characters of the corresponding integral elements are not maximal. 

In the remaining sections of this paper, we shallgive a few singular so- 
lutions for the Riemann-Lanczos problem in 4 dimensions for a number of 
spacetimes and briefly discuss the Riemann-Lanczos problem in 5 dimen- 
sions. 

A. Debever-type Spacetimes: Example 

We expect singular solutions to occur for most of the Debever type space- 
times for which a line element is given by 

ds 2 = dt 2 - f 2 (t, y, z)dx 2 - dy 2 - dz 2 (32) 



with f(t, y, z) being one of the functions given in [12] or in [34]. Specifically, 
for a Debever spacetime with / = y 2 of Petrov type D we are interested in 
singular solutions such that 

£^L a i, c = (33) 

is possible using £^R a bcd = and the fact that V £^ = for £ any Killing 
vector field (=KV). For convenience we shalldenote the local coordinates by 
x 1 := t,x 2 := x,x 3 := y, x 4 := z for this example. Using CLASSI |l]] we 
can determine the dimension of the isometry group which is 4 in this case. 
There are 3 KVs corresponding to the ignorable coordinates t, x and z when 
the above local coordinate frame is used and the components of the fourth 
KV are £ 4 = (z, 0, 0, t). Imposing £^L a \ )C = on the L a b c leads to 



L>tzx 


= 2L>txz 


Ltzy 


= ^Ltyz 


Lfxt 


— Lxzz • 


L t yt 


— Lyzz • 



(34) 

The only non-vanishing independent Riemann-Lanczos equations then are 
AR) AR) AR) . 
Jtxty Jxyxy, Jxzyz- 

Rtxty + ^% y Ltxt — Ptxty = 
— Rxyxy + F xy L xyx — P X yxy = 
Rxzyz Rxzzy xy^xzz • 

One solution which does also satisfy the 6 components of the differential 



) isi 


liven by 






L>txt 


= Ciy 4 , 


Lxyx 


= C 2 y 2 + y 3 


Ltzt 


= c 3 , 


Ltxz 


= C4 , 


Ltyz 


= C 5 y- 2 , 


L X yy 





where C\, ■ ■ ■ , Cq are constants. If we choose C3 = C4 = C% = then only 
3 independent components will remain. On a submanifold with K^L = 
of (^7]) for the above line element this solution corresponds to a singular 
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integral manifold for which the (reduced) characters are sq = s' Q = 3 while 
all higher characters vanish. Such an integral manifold is parameterised 
by the above expressions for the L a b c and their corresponding P a b c d- The 
Vessiot vector fields which span the tangent spaces of this singular solution 
manifold can locally be given by 



V 1 
V 2 
V 3 



V 4 



d_ 

dt ' 

d_ 

dx ' 








dy oLtxt 
d 



2C 2 y- 3 








-— + (2C 3 y + 3y 2 ? 

OLitzy OL, 



dLtyz 

+ 4C iy 3 



4C 22 T 3 




+I2y 



o 



xyx 





UJr txty ur tyzy 







dP f 



+2(C 3 + 3y) 

d_ 

dz ' 







dR 



+ 12y< 







xyxy 



dP x 



ZZl) 



tzyy 



(36) 



where we substituted 



Vi 



txty 



V, 



tyzy 



V: 



xyxy 



6C 2 y" 4 , 
2(C 3 + 3y) , 



(37) 



and where V 3 xty 



T/3 „ nr | t/3 

r xzzy v tzyy 



2V t y zy because of 



). In the next 

section, we are going to dicuss examples of spacetimes where £^L a f, c = 
cannot be implemented but where other singular solutions can occur. 



B. Singular Solution for Godel Spacetime 

Godel spacetime is a perfect fluid spacetime admitting closed time- like curves 
with line element 

ds 2 = a 2 (dt 2 - dx 2 - dz 2 + ^e 2x dy 2 + 2e x dtdy) . (38) 

Again, we are going to replace the local coordinates in such a way that 
x 1 := t,x 2 := x,x 3 := y, x 4 := z. It is well known [25] that this spacetime 



admits a G5 as an isometry group of which 3 Killing vectors C 1 ^ 2 ^ 3 are 
based on the ignorable coordinates t,y, z. The other two are given by 

dx dy' dt dx 2 dy 

But here, the Riemann-Lanczos problem does not admit singular solutions 
which can inherit all the spacetime symmetries. This is the case because 
not all elements of the isometry group of the spacetime are symmetries of 
the Riemann-Lanczos problem. 

However, it is possible to find singular solutions with £^L a b c = imposed 
for the 3 Killing vectors based on ignorable coordinates. An Ansatz where 
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some components are proportional to exponentials e x ,e 2x while all other 
components vanish p3] leads to the solution 



Ltxy 
Ltxt : 



a? e x j a 2 e x 

a e j ^tyx — a e i 



-<xyy 



3a 2 2x 
16 e 



This solution does not satisfy all the differential gauge conditions (||) however 
[ J22| |. This singular solution manifold is again parameterised by the above 
components L a b c and their derivatives P a bcd- Its tangent spaces are spanned 
by 4 Vessiot vector fields which can locally be given as 



V 



V 1 



d_ 

di ' 



d 



8 a 2 

1 e 

dx 8 dL tT 



a 2 
+ —e 







txy 



8 dL 



+ 



3a 2 



+ -e x (— - - - e 2x — 
3 l 16 2 ' " ; 



9 1 . 
e 

16 dPtxyx 3 



e 2 *- 
i 

,5a 
16 







dL 



xyy 
1 

- — e 



2x 



,3a, d 



16 dPtyxx 



V 



v 4 = ^ 



4a 

9z ' 



2x 



9 



an 



xyyx 



(39) 



where it is Vw-r = V 2 xyx . This solution is singular and all Cartan characters 
vanish even sq = s' Q = because no constants remain in the above solution. 

C. Kasner Spacetime 

Kasner spacetime admits a G3 based on the 3 ignorable coordinates x, y, z, 
where we again replaced x 1 := t,x 2 := x,x 3 := y,x 4 := z. If we impose 
£^L a b c = 0, the 6 non-vanishing components of the Riemann-Lanczos equa- 
tions result in 



[I 





— P-txtx 


a 





— P-tyty 


3 





= Rtztz ' 


4 





— Rxyxy 


5 





— Rxzxz 


6 





— P-xzxz 



+ 2Ptxxt 


9 Plr 
* ^ ^txx 


+ %Ptyyt 


? P2 T 
* j. ^tyy 


+ ZPtzzt 


~~ 2—Lt zz 


— 2t 2p2 ~ 


1 P2Ltxx — 



2t 2 ^- x Vx L tyy 
2t 2 ^~ x Vz L txx - 2t 2 ^~ x Vx L tzz 
2t 2 ^- l p 3 L tyy - 2t 2 ^- x p 2 L tzz . 



(40) 



We can now easily see that solving the last 3 equations of fl40|) algebraically 
leaves us with a solution for Lt xx , Ltyy, Lt zz . But, inserting this solution into 
the first 3 equations of (^) and solving for the P a bcd leads to 



Ltxx 
Ltyy 
Ltzz 



■\pi 



\p^-\ 



\P3 



t 2 P3 - 



Ptxxt 
Ptyyt 
Ptzzt 



(I " jPl)Plt 2pi 



lp2)p 2 t 2p2 ' 



h 3 )p3t^- 2 . 
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We see that when we project this down onto the spacetime manifold, we 
obtain that P txxt / L txXtt , P tyy t / L tyyi t and P tzzt ^ L tzZ}t therefore not 
being a solution to (j||). Therefore, for Kasner spacetime no singular solution 
inheriting £^L a i, c = exists. But there are other singular solutions for 
which some of the components L a b c are linear in either x, y or z with no Lie 
symmetries along Killing directions. We decide to make the Ansatz 



Ltxt 


= dt ni x , 


Ltxx 


= C 4 t n * , 


Ltyt 


= C 2 t n *y , 


Ltyy 


= c 5 t" 5 , 


Ltzt 


= C 3 t n *z , 


Ltzz 


= C 6 * n6 , 


Lxyx 


= C 7 t n ?y , 


L>xyy 


= C 10 t nio x 


Lxzx 


= C B t^z , 


Lxzz 


= C u t nil x 


Lyzy 


= C 9 t n *z , 


Lyzz 


= C 12 t n ^y 



where C\,..., C\ 2 and nx,...,ni 2 are arbitrary constants. Inserted into the 
Riemann-Lanczos equations, we obtain 



[T] — Rtxtx — 2Ptxtx + 2Ptxxt — 2^-L txx 



,Pi 
t 

[2] = Rtyty — 2Ptyty + %Ptyyt ~ %^-Ltyy 

[§] = Rtztz — 2Ptztz + 2Ptzzt — 2^-L tzz 

I ^ I = R-xyxy ~ ^Pxyxy ^Pxyyx 2t P2 P 2 Lt xx 2t Pl p\Lfyy 

I 5 I — Rxzxz 2P XZXZ ~\~ 2P XZZX 2t Pi P3^txx 2t ^ x P\Lt zz 

6_ — Ryzyz ^Pyzyz ~i~ 2Py Z zy 2t ^ 3 p 3 Jjtyy 2t ^ 2 p 2 Lf Z z 

= -Rtexy + L tyt t 2pi ~ 1 pi - L xyx t' 1 (p 1 + 2p 2 ) + Pxj/art 

= Rtyxy ~ L txt t 2p2 ~ 1 p2 ~ L xyy r l (2p 1 + p 2 ) + -R^t 

H] = R txxz + L tzt t 2pi - 1 p 1 -L xzx t- 1 (p 1 + 2p 3 ) + P xzxt 

HI] = W 2p3 " 1 P3-ix^ _1 (P3 + 2pi) + P r ^ 

1 = #t raz + W 2P2_1 p2-^ J /^ 1 (P2 + 2p 3 )+ J P^ 

[19] = Rtzyz - Ltytt* 1 *- 1 !* - L yzz t- 1 (p 3 + 2p 2 )+P yzzt , (41) 



where we labelled the equations as explained in Appendix B. Using the above 
Ansatz with n\ = 2p\ — 2,n 2 = 2p 2 — 2,n 3 = 2p 3 — 2, 114 = n\ + 1, n§ = n 2 + 
1, n 6 = n 3 + 1 and n 7 = n 10 = -2p 3 ,n 8 = n u = -2p 2 ,n 9 = n 12 = -2p x , we 
obtain a singular solution for Kasner spacetime with the following constants 
for the above Ansatz 

n. — Cl —El Cr - ° 2 — ei 

° 4 ~~ (Pi-l) 2 ' °5 - (p^x) 2 ' 

n — C3 P3 n — P1C2 

u - 2" ' w - 2-^r ' 

C P1C3 /"( P2C3 



/"I — P2C1 /-I — P3C1 

A computer code, which determines the rather longish expressions for C\,C 2 , 
C 3 in terms of pi,p 2 ,p 3 , can be found in Appendix B or in p^] . For this solu- 
tion all characters and their reduced counterparts vanish identically so that 
so = s' = because all constants are completely determined by pi,p 2 ,p 3 - 
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IV. Comment on the Riemann-Lanczos Problem 
in 5 Dimensions 



On a 5-dimensional spacetime, we can again state the Riemann-Lanczos 
problem. The structure of the Riemann-Lanczos equations remain the same 
as (||) but, we obtain a larger number, namely 50 independent equations due 
to the 50 independent components of the Riemann tensor in 5 dimensions. 
Again, we omit the here 10 components of the differential gauge condition 
L ab s . s = in what follows. 

We also obtain 50 Lanczos components L a b c but imposing the cyclic con- 
ditions (Q) reduces the number to 40 independent components. Looking at 
the dimension of the jet bundle J 1 (M 5 ,R 20 ), we find 5 spacetime dimen- 
sions, 40 independent Lanczos components L abc and 200 independent P a bcd 
totalling N = 245 formal dimensions for the manifold A4. We also find the 
total number of independent 1-forms is given by s < 90, where we have up 

( Ft) 

to 50 independent df^ b( ! d and 40 contact conditions K a b c . The number of 
independent components of a Vessiot vector field V is given by p = N — s 
which here amounts to p > 155. Such a Vessiot vector field can again be 
written as 

V = y£ A + y'Pabce^— + V abcdJ J— (42) 
dx e dL abc dP abcd 

with Vabcd = Vf a bcd}, where {abed} indicates Riemann type symmetries over 
the indices abed. The EDS is again given by (|27]), the only difference being 
the number of equations involved and the range of the indices. 

We wish to find the maximal dimension of a possible involution of Vessiot 
vector fields. The first Vessiot vector field V 1 has at least 155 free ^para- 
metric) components since dim(D) = p > 155. But all additional Vessiot 
vector fields V 2 , • ■ ■ , V 5 need to satisfy a further 40 conditions of the form 

dK abc (V\V j ) = 0, (43) 

where K abc = dL abc —P abce dx e are the contact conditions in local coordinates. 
This leaves the next three Vessiot vector fields with V 2 having at least 
115, V ?J at least 75 and V at least 35 parametric components. The last 
condition for V 5 can cause problems however, as there may not be enough 
free components to form a V 5 in some cases. If we wish to find a suitable 
V 5 , we must have so + si + s 2 + S3 + S4 < 240. We can illustrate this as 
follows: 



s < 90 


N 


— sq > 155 , 


s + si < 130 


N 


— so — si > 115 , 


so + si + s 2 < 170 


=> N 


- so — si — s 2 > 75 , 


so + si + s 2 + s 3 < 210 


=> N 


- s + si + s 2 + s 3 > 35 , 


so + si + s 2 + s 3 + s 4 < 250 


=> N 


- so — si - s 2 — S3 - S4 > -5 



where the second column gives the number of free components for V 1 ^ 2 
and so on. Note that the fact that the sum of all free components 155 + 115 + 
75 + 35 = 380 > 245 is permissible because 380 is the minimal total number 
of free but not necessarily independent components for all V , • • • , V 5 . 

If we abandoned the cyclic conditions for a moment, the numbers would 
be: N = 250, s < 100, p > 150. Then V 1 would have at least 150 components 
left, V 2 at least 100 components, V 3 at least 50 components and V 4 would 
have no parametric components left in the worst case when all the Si are 
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maximal. There would not be any parametric components left to form V 5 . 
Therefore, in the generic case no 5-dimensional involution of Vessiot vector 
fields exists. This means that in the generic case we cannot even form a 
singular 5-dimensional chain of integral elements but only a 4-dimensional 
at most. Therefore, neither the first nor the second of the Cartan-Kahler 
theorems hold in this case. 

Apart from that, we expect 20 identities of the same kind as the identities 
in H to occur. Using the computer code in Appendix C given in p|| adapted 
to the 5-dimensional case, we could find the following results for the reduced 
Cartan characters for a number of spacetimes such as Minkowski space, 
the 5-dimensional conformally flat spacetime and the Debever spacetimes. 
Their reduced characters (s^, s' 2 , s' 3 , S4, s' 5 ) when not including the differential 
gauge condition (||) were (40,39,35,26,10) and the second set including the 
differential gauge conditions (Q) was (40,39,35,26,0). The computer code 
also confirmed that neither of the systems were in involution. We conclude 
that in order to be a viable problem in 5 dimensions, the Riemann-Lanczos 
problem needs to be modified significantly. 



Conclusion 

We gave a review of the theory of exterior differential systems and Pfaffian 
systems in particular. We illustrated this theory with the 2-dimensional 
example of a coordinate transformation expressed as a complete Pfaffian 
system. 

Singular solutions for the Riemann-Lanczos problem in 4 dimensions do 
exist and we found some for Kasner, Godel and for a Debever type spacetime. 

For the Riemann-Lanczos problem in 5 dimensions, 20 internal identities 
of the same type as those in 4 dimensions can occur as a computer code 
suggests and a 5-dimensional involution of Vessiot vector fields cannot be 
achieved in a generic case. 
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Appendix A: The Cartan-Kahler Existence Theo- 
rems 



The Cartan-Kahler existence theorems [|], 45, 27, 10] are stated here. The 
first Cartan-Kahler theorem specifies under which conditions a unique inte- 
gral manifold of dimension {p + 1) can be found from a p-dimensional one: 



Theorem 7 First Cartan-Kahler Existence Theorem 

Let M be a formal M- dimensional manifold andM be a p- dimensional inte- 
gral manifold with a regular integral element T x (AT) , where (E p ) x = T X (N) at 
a point x on M . Further, there exists a submanifold T , where N C T C Ai 
such that dim^) = M — i p +i, dim^^T) n H({E P ) X )) = p + 1, where 
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tp+i = dim(H((E p ) x )) — p — 1. Then, there exists a unique integral manifold 
X around x such that dim(X) = p + 1 and J 7 D X D N . 

This means that for a given p + 1-dimensional manifold J-, we are looking 
for the submanifold X C J- such that dim(A') = p + 1 and the theorem 
tells us that X exists and is unique. Here, the integer t p+ \ is defined as 
t p+ i = dim(E p ) x —p—1, < p < n and s p = t p — t p+ \ — 1. For to we have 
to = Y^i=Q s i + n = dimTiT^o) . For the second Cartan-Kahler existence theo- 
rem, we need our EDS E to be given in normal form which in this context 
means that there exists a local coordinate system (x 1 , .., x n , y 1 , .., y r ), where 
r = N — n, such that: 

i) J(S ) is defined by y*o-n+i = . . . = y iV-n = j 

ii) H{{E p ) x ) = {d x i,-- ■ ,d x n,d y(so+ ... +Sp _ 1+ i), ■ ■ ■ ,d y ( t0 -n))}, 0<p<n, 
Hi) (E p ) x = (d x i, ■ ■ ■ ,d xP ), 1 < p < n , 

iv) The integral point (E°) x is given by(E°) x = (0, • • • , 0) , so that the point 
x coincides with(0, • • • , 0)(iV-times) on M. . 

Then, we can state the second Cartan-Kahler existence theorem: 

Theorem 8 Second Cartan-Kahler Existence Theorem 

Given a regular chain of integral elements (E°) x C • • • C (E n ) x of an EDS 

E given in normal form. Consider a set of initial data 

/l) • • " j /s ; ( s arbitrary constants), 

f S()+ i(x 1 ), • • • , / S0+S1 (x 1 ) , (si arbitrary functions of 1 variable each), 

fs +-+s n -i+i{x 1 ,---,x n ),---,f to -n{x 1 ,---,x n ) , (s n arbitrary functions of 
n variables x 1 , ■ ■ ■ , x n ) . 

For sufficiently small values of all the above fi and their first-order deriva- 
tives there exists a unique integral manifold defined by y l (x 1 , ■ ■ ■ , x n ), yi = 
0,1 <i < t -n<j < N-n, such that y^x 1 , ■ ■ ■ , x p , 0, ■ ■ • , 0) = f^x 1 , ■ ■ ■ , x p ) 
for so + • • • s p-i < i < s o + • • • + s p , < p < n. 

The highest non- vanishing Cartan character say Sk, k < n is an invari- 
ant of an EDS and its value gives the number of arbitrary functions of k 
variables involved in the general solution. The other characters have the 
same meaning only if we used normal local coordinates in the above 
sense. 



Appendix B: Labelling of the f^ d 
for Kasner Spacetime 



and Constants 



The labelling of the equations for the Riemann-Lanczos problem in 4 dimen- 
sions as used in section III.C was carried out in the following way: 



JR) 
J1212 > 



JR) 



^ Jl313 ' 



AR) 
/1414 i 



AR) 
/2323 ' 



(R) 



13 



2424 ' 



AR) 

J 1323 ' 

_ AR) 

^ J 1224 ' 



10 



14 



AR) 
'3434 > 

AR) 

J 1324 ' 

AR) 

J 2324 ' 



(R) 
1213 ' 



HH ^ /l42 I 



AR) 
J 1214 > 



e(R) 



12 



16 



AR) 

J 1223 ' 
^ /1314 ) 
^1234 ' 



17 



f (R) 

^1334 ' 



18 



J 2334 ' 



19 



AR) 

^1434 > 
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AR) 

J 2434 • 
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Next, we would like to give the computer code used to determine the con- 
stants for the singular solution for Kasner spacetime. The expressions for 
the constants C4, ■ • • , C\i given in section III.C above can also be found in 
||22|| . They will be used as input for the REDUCE code which determines 
Ci, C2, C3 completely in terms of pi,P2,P3- This code is given by: 

°/ Constants for a singular solution for the Riemann-Lanczos 
°/ problem for Kasner spacetime: 
°/ pl+p2+p3=l ; pl**2+p2**2+p3**2=l ; 

/oFirst we solve equations (1) to (3) and (7) to (12) and obtain: 

C4:=Cl/(pl-l)-pl/2; 

C5:=C2/(p2-l)-p2/2; 

C6:=C3/(p3-l)-p3/2; 

C7:=(C2*pl/(2*p3+2*p2+pl)) ; 

C8:=(C3*pl/(2*p3+2*p2+pl)) ; 

C9 : = (C3*p2/ (2*p3+2*pl+p2) ) ; 

C10:=(-Cl*p2/(2*pl+2*p3+p2)) ; 

Cll:=(-Cl*p3/(p3+2*p2+2*pl)) ; 

C12:=(-C2*p3/(p3+2*p2+2*pl)) ; 

°/ We write (4), (5) and (6) for the remaining CI, C2 and C3 and 
°/ solve : 

s4 : =2*C7-2*C10+2*p2*C4+2*pl*C5+pl*p2 ; 

s5:=2*C8-2*Cll+2*p3*C4+2*pl*C6+pl*p3; 

s6 : =2*C9-2*C12+2*p3*C5+2*p2*C6+p2*p3 ; 

solve (s4=0, CI) ; 

Cl:=rhs first ws; 

solve (s5=0,C2) ; 

C2:=rhs first ws; 

solve(s6=0,C3) ; 

C3:=rhs first ws; 

off nat ,nero , echo ; 
out "H : /constants" ; 
C1:=C1; C2:=C2; C3:=C3; 
shut "H : /constants" ; 
on nat ,nero , echo ; 
end; 

The computer code above produces an output file called "constants" with 
the constants C±, C2 and C3 expressed in terms of the pt. 



The term torsion in the context of EDS has nothing to do with the torsion 
occurring in the theory of affine connections. 
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